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Abstract It has been suggested that the frequency in the 
co-rotating innermost stable circular orbit (ISCO) about 
a compact stellar remnant can be determined through X- 
ray observations of low-mass X-ray binaries, and that its 
value can be used to constrain the equation of state of 
ultradense matter. Upon constructing numerical models 
of rapidly rotating strange (quark) stars in general rel- 
ativity, we find that for stars rotating at the equatorial 
mass-shedding limit, the ISCO is indeed above the stellar 
surface, for a wide range of central energy densities at a 
height equal to 11% of the circumferential stellar radius, 
which scales inversely with the square root of the energy 
density, po, of self-bound quark matter at zero presure. 
For these models, the ISCO frequency is 81.5±1.5% of the 
stellar rotational frequency, whose maximum value fix = 
\/3.234 GpQ is attained for a model close to the maximum- 
mass model, with Af = 2.86M0(/9o/4.2xlOi'*gcm-3)-i/2 

In contrast to static models, ISCO frequencies below 1.1 
kHz are allowed — in fact, at the canonical value po ~ 
4.2 X lO^'^gcm"'^, the ISCO frequencies of rapidly rotat- 



ble orbit about the star ( 


Kluzniak and Wagoner, 1985 


pyunyaev and Shakura, 1986; 


Kluzniak et al., 1990) leads 



to significant constraints on the equation of state of mat- 
ter at such densities, as very few models of ultra-dense 
matter admit neutron stars of mass high enough to al- 
low maximum orbital frequencies as low as the observed 
values in the quasi-periodic oscillations (QPOs) — for in- 
stance, th e QPO frequency in 4U 1820-30 saturates at 



1.07 kHz (IZhang et al., 1998|) 



ing strange stars can be as low as 0.9 kHz for a 1.3Mq 
strange star. Hence, the presence of strange stars in low- 
mass X-ray binaries cannot be excluded on the basis of 
the currently observed frequencies of kHz QPOs, such as 
the cut-off frequency of 1066 Hz in 4U 1820-30. 



Key words: dense matter - equation of state - stars: 
binaries: general - X-rays: stars 

1. Introduction 

The discovery of millisecond variability in the flux of low- 
mass X-ray binaries (LMXBs) has raised the prospect of 
constraining the properties of matter at supranuclear den- 
sities, which is thought to make up the compact stellar 
remnant in these sources. Conventionally, the compact ob- 
jects is taken to be a neutron star, and it has been shown 



Still, the evolutionary status of LMXBs and the nature 
of the accreting compact object are not clear. It is known 
that the X-ray bursters cannot be black holes, because 
their photospheric radius and the temperature during the 
burst both tend to a definite value, thus showing the pres- 
ence of a stellar surface, which is also required to explain 
the (type I) X-ray bursts as thermonuclear explosions of 
accreted material. The inferred radii (and masses) are con- 
sistent with models of neutron stars, but it is possib le that 
the compact o bject is a "strange" , i.e. quark, star ( [Cheng 
and Dai, 1996| ). If it were, at least in the sources 4U 1820- 



30 and 4U 1636-53, then the energy density of self-bound 
quark matter at zero pressure would have to have the un- 
usually low value po < 4.2 x lO^'^gcm"^, if the maximum 
observed frequencies of the kHz QPOs were the orbital fre- 
quencies in the ISCO about slowly rotating strange stars; 
and the observed QPO frequencies could not be the or- 
bital frequencies at the surface of such stars, for any value 
of Po dBulik et al., 1999| ). 

However, it seems likely that in these very old LMXBs, 
the compact star has been spun up through accretion 
to very highQ frequencies — a neutron star would have 
to accrete only ~ 0.2Mq to attain angular momentum 
J = O.GGAP/c ( [Kluzniak and Wa goner, 1985[ ). A strong 
magnetic field could prevent rapid rotation, but strange 
stars may not support such a field (Horvath 1999). We 



(Kaaret et al., 1997; Kluzniak, 1998) how the assumption 



that the highest observed frequency in the X-ray flux is 
the orbital frequency in the innermost (marginally) sta- 



^ The frequency of ~ 300 Hz observed during some X-ray 
bursts is unlikely to be the rotational frequency of the star, 
as the long duration of the oscillation (~ 10 s), its harmonic 
content, and its variation during the burst, seem incompatible 
with the flux modulation being caused by a "hot spot" on the 
stellar surface. 
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Figure 1. Radius vs. mass for strange stars (near their 
maximum mass). Models of both static (thin hnes) and 
maximally rotating (thick lines) stars are shown for three 
values of the energy density of quark matter at zero 
pressure, po- 4.2 x lO^^gcm^'^ (continuous line) 5.3 x 
lO^^'gcm"^ (long-dashed line), 6.5 x lO^^gcm"'^ (dashed 
line). The maximum-mass models for these values of po 
are indicated by a filled circle, an open circle and a star, 

^1/2 



Figure 2. Mass vs. central energy density for the same 
models as in Fig. 1. For each sequence, a star is stable 
only at values of p smaller than, approximately, the one 
corresponding to the maximum mass. 



-1/2 



and is 2Mq for 



4x10^^ gem ^ . D etailed models of str ange st ars have 



respecii ivcly. Note that the radiua and maaa acalc aa 



A larg e incrcaae of the radiua and maximum maaa ia cvi — 
dent as the stellar rotation rate increases from zero to the 
equatorial mass-shedding limit. 

show that if LMXBs harbour rapidly rotating strange 
stars, the constraints from kHz QPOs on the stellar mass 
and on po are relaxed to a surprising degree. We also dis- 
cuss the rotational frequency of "Keplerian" models of 
strange stars, i.e. ones in which the rotation rate at the 
stellar equator is equal to the orbital frequency for the 
star, at the same equatorial radius. 

2. Strange stars 

Quark stars are likely to exist if the ground state of mat- 



of a static strange star scales as Pq 

been constr ucted by ( Alcock et al., 1986 ) and ( Haense] 
et aL, 1986| ). 



Astrophysical implications of these ideas are not yet 
clear. It has been pointed out that young, glitc hing radio 
pulsars are probably neutron stars ( Alpar, 1987 ), and that 
strange stars are unlikely to be present in Hulse- Taylor 
type binaries, as their coalescence may lead to dispersal of 
nuggets of quark matter which would have precluded the 
formation of young neutron stars in the Galaxy (Madsen 
1988, Caldwell and Friedman 1991). However, there seems 
to be no objection to millisecond pulsars or the compact 
objects in LMXBs be ing strange stars ( Kluzniak, 1994 ; 
[Cheng and Dai, 1996| ), and it has even been suggested 
(Madsen 1999) that millisecond pulsars can be formed di- 
rectly in supernovae, if they are strange stars [unlike neu- 
tron stars, whose rotation rate would be quickly damped 
by the r-mode instability ( Andersson et al., 1999D 



ter at l arge atomic number is in the form of a quark fluid, 
which would then necessarily be composed of about equal 



bom et. al., 199q ), ( Kokkotas and Stergioulas, 1999| )] 



indl- 



numbers of up, down, and strange quarks (Bodmer, 1971). 
Today, such matter is called strange matter. Its thermo- 
dynamic properties have been discussed in detail by Farhi 
and Jaffe (1984) in the context of quantum chromody- 
namics. The first relativistic model of stars composed of 



Our work is informed by the question whether the pres- 
ence of strange stars in LMXBs can be excluded on the 
basis of the observed timing properties of these sources 
Specifically, it has been asked whether the observed fre- 



quark matter was computed by (Brecher and Caporaso, 



1976). The cosmological consequences of the presumed ex- 
istence of strange matter were first discussed in detail by 
( Wittcn, 1984 ), who also showed that the maximum mass 



^ In the conversion of accreted matter to quark matter, sub- 
stantial energy will be released at the base of the crust. How- 
ever, since the mass accretion rate in LMXBs is only inferred 
from the X-ray luminosity, there is no direct way of differenti- 
ating this from the gravitational binding energy release, if the 
conversion of matter is proceeding at a quasi-steady rate. 
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Table 1. Keplerian models of strange stars, for po = 4.2 x lO^^'gcm 





M 


R 


Q. 


n+ 


h+ 


h- 


J 


I 




T/W 




Zb 




lO'^gcm-3 




km 


lO^'s-i 


lO^s-i 


km 


km 


GMl/c 


10*'^ gem 


2 










5.000 


1.336 


16.81 


7.049 


5.664 


1.86 


9.32 


2.349 


2.929 


.3192 


.2785 


.2731 


.8931 


.2821 


5.250 


1.598 


17.29 


7.257 


5.877 


1.90 


11.2 


3.123 


3.782 


.3365 


.2709 


.3279 


1.045 


.2942 


5.500 


1.800 


17.62 


7.432 


6.051 


1.87 


12.7 


3.760 


4.447 


.3500 


.2645 


.3739 


1.176 


.3048 


5.750 


1.963 


17.75 


7.591 


6.217 


1.86 


13.9 


4.299 


4.978 


.3630 


.2590 


.4145 


1.294 


.3113 


6.105 


2.161 


17.89 


7.786 


6.389 


1.88 


15.6 


5.003 


5.647 


.3750 


.2534 


.4676 


1.453 


.3199 


6.500 


2.304 


17.87 


7.979 


6.582 


1.84 


16.7 


5.474 


6.030 


.3906 


.2465 


.5121 


1.589 


.3260 


6.600 


2.336 


17.87 


8.012 


6.615 


1.84 


16.9 


5.590 


6.131 


.3930 


.2457 


.5221 


1.619 


.3270 


6.750 


2.388 


17.86 


8.080 


6.678 


1.84 


17.4 


5.788 


6.296 


.3969 


.2441 


.5393 


1.674 


.3291 


7.275 


2.520 


17.79 


8.280 


6.864 


1.82 


18.6 


6.258 


6.643 


.4094 


.2387 


.5875 


1.829 


.3348 


8.668 


2.730 


17.51 


8.708 


7.229 


1.79 


20.5 


7.007 


7.072 


.4312 


.2285 


.6835 


2.153 


.3458 


9.023 


2.753 


17.40 


8.803 


7.327 


1.76 


20.7 


7.055 


7.043 


.4367 


.2257 


.6999 


2.207 


.3474 


9.658 


2.796 


17.23 


8.958 


7.459 


1.75 


21.1 


7.180 


7.045 


.4437 


.2221 


.7290 


2.309 


.3501 


10.33 


2.823 


17.06 


9.103 


7.582 


1.74 


21.3 


7.233 


6.983 


.4504 


.2186 


.7535 


2.394 


.3523 


11.06 


2.844 


16.89 


9.245 


7.699 


1.73 


21.5 


7.259 


6.901 


.4562 


.2154 


.7760 


2.474 


.3544 


11.84 


2.854 


16.69 


9.385 


7.823 


1.72 


21.6 


7.232 


6.772 


.4625 


.2119 


.7953 


2.541 


.3555 


12.67 


2.861 


16.52 


9.522 


7.933 


1.71 


21.7 


7.211 


6.656 


.4672 


.2090 


.8140 


2.608 


.3573 


13.56 


2.860 


16.33 


9.648 


8.047 


1.69 


21.7 


7.124 


6.489 


.4727 


.2056 


.8271 


2.652 


.3579 


14.51 


2.852 


16.14 


9.782 


8.153 


1.69 


21.7 


7.038 


6.323 


.4773 


.2026 


.8415 


2.702 


.3589 


15.53 


2.842 


15.97 


9.909 


8.256 


1.68 


21.6 


6.937 


6.153 


.4812 


.1998 


.8532 


2.741 


.3600 


16.63 


2.829 


15.78 


10.03 


8.363 


1.67 


21.5 


6.829 


5.980 


.4852 


.1970 


.8638 


2.777 


.3606 



quency of the kHz QPOs in sources such as 4U 1820 -30 
(1.07 kHz) is not too low to be compatible with the max- 
imum mass of strange stars (Bulik et al., 1999). For this 
reason, in constructing our models of strange stars, we 
have focussed on an equation of state which yields the 
largest masses of static strange star models: 



P = (p-po)/3. 



(1) 



We have also investigated the more general case, where the 
factor of 1/3 in eq. (1) is replaced by a different positive 
constant, a <1. 

3. Keplerian models of strange stars 

We have computed exact numerical models of strange 
stars in general relativity using the (Stergioulas and Fried- 
man, 199^ ) code (see Stergioulas 1998 for a description) 
In this code, the equilibrium models are obtained follow- 
ing the KEH (Komatsu et al. 1989) method, in which the 
field equations are converted to integral equations using 
appropriate Green's functions. 

The detailed expected properties of strange stars de- 
pend on the adopted theory of interactions. All models 
presented here were constructed using eq. (1) for the equa- 
tion of state. For our models of rotating stars, we find that 
mass and radial quantities (e.g., the stellar radius and the 

— 1/2 

height of the ISCO above it) scale as Pq , just as for the 

1 /2 

static stars, while the frequencies scale as . For the 

more general e.o.s. P = a(p-po), we con firm the approx- in rapidly rotating stars (^hibata and Sasaki, 199^). 



as a^/**, and the rotation rate of the star scales as a~^/*. 
Thus, between the scalings with po and with a, the nu- 
merical results presented here can at once be extended to 
the general e.o.s. of strange matter. 

In Figs. 1 through 3, we present the mass, radius and 
the ISCO frequencies in our Keplerian models, for three 
values of and compare them with the static models. 
In Fig. 4 we present the ISCO angular frequencies as 
a function of the central energy density of the strange 
star, and also exhibit the (larger) angular frequency of 
the star itself. Note, that because of the scalings with en- 
ergy density, the maximal rotation rate of a strange star 
(which is very close to the rotation rate of the maximum- 
mass model in Table 1) is described by the simple formula 
Q,K — -\/3-234 Gpo: where G is the gravitational constant. 

Table 1 presents in detail, the various stellar param- 
eters obtained in our calculation for Keplerian strange 
stars, modeled with eq. (1) for the value po = 4.2 x 
lO^^gcm"'^. Note that in all cases, the co-rotating ISCO 
is above the stellar surface and, as exhibited in Fig. 3, the 
ISCO frequencies are much lower for the Keplerian mod- 
els than for static models (at fixed stellar mass) and differ 
considerably from their lowest-order slow-rotation approx- 
imation. The significant departure from the slow-rotation 
result is explained by the unusually large oblateness of 
rapidly rotating strange stars, and the fact that the ISCO 
frequency and height depend not only on the angular mo- 
mentum, but also on the stationary quadrupole moment. 



imate scalings of (Lattimcr et al., 1990) to within 9 %, 



that is, we find that for our Keplerian models, the max- 
imum stellar mass scales as a}/"^, the stelar radius scales 



4 



Stergioulas, Kluzniak, Bulik: Keplerian frequencies of strange stars 



1.6 



1.2 



- 

1 


- 

1 1 
\ 


- 


\ 

\ 




\ 




\ 

^ \ 




! \ 








\ 




\ 




f \ - 

/ \ 




/ \ 




\ ^ ^ J 


1 , 





1.5 



2.5 



1.4 



1.2 - 



1 - 



0.8 - 




101' 



M[M^ 



2x1015 
p[ gcm-3] 



3x1015 



4x1015 



Figure 3. The frequency of the co- rotating innermost sta- 
ble circular orbit as a function of mass for static models 
(thin, continuous line) and for strange stars rotating at the 
equatorial mass-shedding limit (thick lines, in the style of 
Fig. 1). For the static models, /+ = 2198 Hz(M0/M), and 
the minimum ISCO frequency corresponds to the maxi- 
mum mass, denoted by a filled circle, an empty circle, and 
a star, respectively for /9o/(10^'*5 cto~^) = 4.2, 5.3, and 
6.5. Note that the ISCO frequencies for rapidly rotating 
strange stars can have much lower values, and /_)_ < 1 kHz 
can be achieved for strange stars of fairly modest mass, 
e.g. 1.4M0, if the star rotates close to the equatorial mass- 
shedding limit. 



4. Conclusions 

We have calculated exact models of rotating strange 
stars — these computationss are in excellent agreement 
(Table 2) with the very recent results obtain ed by a highly 
accura te co de based on spectral methods ( Gourgoulhon 
L999| ). We found the scahngs of A/, R, n with the 



et al. 

parameters a, and po in the equation of state of self-bound 
quark matter P = a{p — po). In addition, we calculate 
the innermost stable orbits and find that, unlike for static 
models, for strange stars rotating at the equatorial mass- 
shedding limit the orbital frequencies can extend to values 
below 1.1 kHz. For the same models the radius of the ISCO 
is about 11% larger than the circumferential stellar radius, 
independently of the central density, or the value of po- 

Our results show that the highest observed QPO fre- 
quencies in low-mass X-ray binaries (such as 1.07 kHz in 
4U 1820-30) could be the orbital frequencies in the in- 
nermost stable circular orbit about strange stars, if only 
the stars are rotating sufficiently rapidly (as is expected 
in these old accreting sources). Thus, the compact objects 



Figure 4. The angular frequency of strange stars (thick 
lines) rotating at the equatorial mass-shedding limit, and 
of their co- rotating ISCOs (thin lines), as a function of 
central energy density. The style of lines corresponds to 
the same values of po as in Fig. 1. 

Table 2. Comparison of two codes. 

This work Gourgoulhon diff. 











et al. 


[%] 


Pc 




cm ^ 


1.261 


1.261 




M 


Mq 




2.83392 


2.831 


0.1 


R 


km 




16.4252 


16.54 


0.7 


n 


lO^s" 


1 


9.56261 


9.547 


0.2 






1 


9.57637 


9.547 


0.3 


T/W 






0.210091 


0.210 


0.04 


cJ/GMl 






7.09331 


7.084 


0.1 


I 


10«g 


cm^ 


6.51906 


6.534 


0.2 


Zp 






0.807995 


0.8070 


0.1 


Zb 






2.58820 


2.584 


0.2 








0.356272 


0.3608 


1.2 








0.466000 


0.4618 


0.9 



in LMXBs could, in principle, be rapidly rotating strange 
stars. Further conclusions about the nature of LMXBs and 
of the kHz QPOs would be possible, if either the mass or 
the rotational period of the accreting stellar remnant were 
known. 



This research was supported in part by the NSF grant 
PHY94-07194. 
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